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ABSTRACT 

Context. We present VLT-ISAAC medium resolution spectroscopy of the HH34 and HH1 jets. 

Aims. Our aim is to derive the kinematics and the physical parameters and to study how they vary with jet velocity. 

Methods. We use several diagnostic lines, such as [Fell] 1.644/jra, 1.600yum and H 2 2.122yum, in order to probe the atomic and 

molecular components. 

Results. In the inner jet region of HH34 we find that the atomic and molecular gas present two components at high and low velocity 
(the so-called HVC and LVC). The [Fen] LVC in HH34 is detected up to large distances from the source (> 1000 AU), at variance 
with TTauri jets. In Hi 2. 122 pm, the LVC and HVC are spatially separated, with an abrupt transition from low to high velocity 
emission at ~ l'.'S. We moreover detect, for the first time, the fainter red-shifted counterpart down to the central source. In HH1, we 
trace the jet down to ~ 1" from the VLA1 driving source: the kinematics of this inner region is again characterised by the presence 
of two velocity components, one blue-shifted and one red-shifted with respect to the source LSR velocity. We interpret this double 
component as arising from the interaction of two different jets. We suggest that the red-shifted component could be part of the HH501 
jet. Electron densities and mass fluxes have been measured separately for the different velocity components in the HH34 and HH1 
jets. In the inner HH34 jet region, n e increases with decreasing velocity, with an average value of ~ 1 x 10 4 cm~ 3 in the HVC and 
~ 2.2 x 10 4 cirT 3 in the LVC. Up to ~ 10" from the driving source, and along the whole HH1 jet an opposite behaviour is observed 
instead, with n e increasing with velocity. In both jets the mass flux is carried mainly by the high-velocity gas: lower limits on the mass 
flux of 3 - 8 x 10~ 8 M yr~' have been found from the luminosity of the [Fen] 1. 644 pmi line. A comparison between the position 
velocity diagrams and derived electron densities with models for MHD jet launching mechanisms has been performed for HH34. 
While the kinematical characteristics of the line emission at the jet base can be, at least qualitatively, reproduced by both X- winds and 
disc-wind models, none of these models can explain the extent of the LVC and the dependence of electron density with velocity that 
we observe. It is possible that the LVC in HH34 represents denser ambient gas entrained by the high velocity collimated jet. 

Key words, stars: circumstellar matter - Infrared: ISM - ISM: Herbig-Haro objects - ISM: jets and outflows - ISM:individual: 
HH34, HH1 



1. Introduction 

Protostellar jets and accretion discs are phenomena intimately 
connected. The relation between both structures and the jet laun- 
ching mechanism is, however, still not well understood. Several 
models have been developed to constrain the physical mecha- 



nism by which mass is accelerated from the star's vicinity ( |Shu| 



|et al. 


1995 Ferreira 1997|Matt & Pudritz 


2005) by the ac- 


tion of so-called "rotating magnetospheres" 


Camenzind 1990). 



Observations of the inner part of the jet structure are required, 
however, in order to verify the predictions of such kinds of mod- 
els. In particular, NIR spectroscopy is an important tool to inves- 
tigate the jet structure closest to the driving source. This is espe- 
cially true for jets in Class 0/1 objects, where the high extinction 
prevents us from observing the regions close to the source with 
standard optical tracers. Several observational studies have in- 
deed been carried out during recent years employing near-IR line 



Send offprint requests to: R. Garcia Lopez, e-mail: 
garcia@mporzio . astro . it 

* Based on observations collected at the European Southern 
Observatory, La Silla, Chile (ESO programmes 0.74.C-0286(A)). 



diagnostics on jet beams from young embedded sources (e.g. 
Nisin i et al.|2005t |Davis et al.|2006l |Podio et al.|2006l |Takami| 
|et al.|2006| ). In these studies the main lines investigated are [Fe n] 
lines (e.g. 1.644 yt/m and 1.600/im) and H2 lines (the 2.12//m 
line is among the brightest lines). These works have addressed 
both the physical parameters and the kinematical properties of 
the sources, and their similarities with more evolved Classical 
T Tauri Stars (CTTSs). Kinematic ally, both Class I and CTT 
stars, show a similar behaviour, with the presence of a colli- 
mated and large scale jet at high velocity (~200-400kms~ I ), 
called the High Velocity Component (HVC), and a compact gas 
component localised around the central source at lower velocity 
(~0-50kms~' ) named Low Velocity Component (LVC). This 
supports the assumption that the accretion and ejection in Class 
I objects occurs by the same mechanism as for T Tauri stars. 
Synthetic Position Velocity (PV) diagrams constructed for diffe- 
rent classes of MHD jet models (i.e. disc-winds, |Garcia et al. 
|2001| and X-winds, Sha ng et al.|2002| ) predict the presence of di- 
fferent velocity components at the jet base, although the details 
of the observed kinematical features are still not well reproduced 
by these models. 
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Table 1. List of detected lines 



HH34 short K 











Line id. 






Jet 



"7/2 

[Fen]a 4 Z>7/2 - a 4 F(>/2 
[Fen]a 2 P V2 -a 4 P 3/ 2 
[Fen]a 2 H n /2 - a 2 G 9/2 
[Tin]a 2 F 5/2 - a 4 F 3/2 
[Tiii]a 2 F 7/2 - a 4 F 5/2 
H 2 1-0 S(l) 
H 2 2-lS(2) 
H 2 1-0 S(0) 
H 2 2-l S(l) 
Bry 



1.6440 
2.1333 
2.2244 
2.1605 
2.0818 
2.1218 
2.1542 
2.2235 
2.2477 
2.1661 



HH34, HH1 
HH34, HH1 

HH34 

HH1 

HH1 
HH34, HH1 

HH1 

HH34 

HH34 

HH1 



Vacuum wavelengths in microns. 



On their physical properties, however, jets from Class I and 
CTTSs show some differences: in general, densities and mass 
fluxes derived for jets from Class I sources are higher than those 
in jets from CTTSs, as one would expect in sources with higher 
accretion rates 



embedded in a dense environment (e.g. Davis 
|et al.||20Q3| |Antoniucci et al.||20"0"8) . Moreover, the molecular 
component traced by H 2 NIR lines is significantly enhanced in 
the younger sources, and in some of them the jet shows up only 
in molecular form. 

In this framework, an important observational test is to mea- 
sure the relevant physical parameters separately in the different 
velocity components of the flows. This kind of diagnostics can 
show if the differences in the physical properties between jets 
from Class I sources and CTTSs persist in both the HVC and 
LVC, and provides further constraints on the origin of these lat- 
ter. Differences in the excitation conditions of the various velo- 
city components have been measured in jets from CTTSs. Such 
a kind of study has, however, never been attempted for jets from 
Class 0/1 stars. 

Here, we present recent results on the HH34 and HH1 jets, 
two classical outflows that are bright in the near-IR (e.g. [Davis] 



etal. 2000 2001 



Eisloffel et al. 2000 1. Their driving sources are 



HH34 IRS and VLA1, two young embedded Class 0/1 objects 
situated at a distance of 450 pc and 460 pc, respectively, in a high 
phase of accretion (Antoniucci et al.||2008| |Chini et aL 2001 



Velocity resolved, near-IR spectroscopy of HH34, addressing 
the kinematics of the inner region, has been presented by |Davis| 
[etaE1 ( r200T1 ), |Davis et al.| ( [2003] > and |Takami etd^ l |2006^ . More 
recently, we have observed the HH34 and HH1 jets in low reso- 
lution spectroscopy from 0.6 to 2.4 /vm, and derived the relevant 
physical parameters of the jets, as a function of the distance from 
the exciting source, adopting an analysis combining optical and 
near-IR line ratios ( |Nisini et al.|2005[|Podio et al.|2006) . In this 
work, we present velocity-resolved NIR spectra obtained with 
the spectrograph ISAAC on the Very Large Telescope (VLT) of 
the European Southern Observatory (ESO). This instrument co- 
vers important diagnostic lines of the ionised ([Fe n] 1.644, 1.600 
yum) and molecular (H 2 2. 122 /im) gas. Our aim was studying the 
kinematics of the atomic and molecular components, and to de- 
rive the electron density and mass ejection rates in the different 
velocity components. The paper is structured as follows: the ob- 
servations and data reduction methods are described in Section 
2. In Section 3, we present the results on the kinematics of the 
[Fen] and H 2 lines. In Section 4, we derive the electron density 
and mass flux from the [Fe n] line luminosity and ratios, and dis- 
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Fig. 1. Continuum-subtracted K-band spectral images of the 
HH34jet. 



cuss the variation of these parameters as a function of velocity 
and distance from the central object. In Section 5, we discuss the 
comparison of our results with the predictions of proposed jet 
launching models. Finally, we draw our conclusions in Section 
6. 



2. Observations and data reduction 

Observations were obtained on 28-29 December 2004 at the 
ESO VLT telescope on Cerro Paranal, Chile, using the infrared 
spectrograph and camera ISAAC at medium resolution. The 
adopted slit width is 0/3, corresponding to a nominal resolu- 
tion R-8900 in the K-band and 10000 in the H-band. The spa- 
tial scale of the camera is O'.'146/pixel. We have taken H (1.57- 
1.65 //m) and K (2.07-2.20 //m) band spectral segments with the 
slit aligned along the jets. An additional segment in K-band, 
covering the range 2.19-2.31 /mi, was also acquired on HH34. 
The adopted position angles (PA) were 145° and -15° for the 
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Fig. 2. Continuum-subtracted K-band spectral image of the HH1 jet. 



2.16 



HH1 and HH34 jets, respectively. Total integration times were 
3300 s in the K-band for both the jets and 5100 s and 6900 s 
in the H-band for HH1 and HH34, respectively. Data reduc- 
tion was performed using standard IRAfQ tasks. The IR spec- 
tra were corrected for the atmospheric spectral response by di- 
viding the object spectra by the spectrum of a B5III star show- 
ing telluric lines. Wavelength calibration was performed using 
the atmospheric OH emission lines. The H and K-band effec- 
tive resolution measured from the OH lines is ~8600 and ~7800 
for both objects. Flux calibration was performed using a photo- 
metric standard star. The standard stars have been observed at 
similar air mass and seeing conditions as the objects, therefore 
no corrections for flux losses due to the slit width smaller than 
the seeing (~(X'6 in both the nights) have been applied. This pro- 
cedure probably leads to some flux losses in the more external 
section of the jet (> 15" from the source) where the intrinsic jet 
diameter, as measured from HST images (Reipurth et al.|[2002 



widens to more than 076. In the case of HH34, we subtracted the 
continuum emission from the HH34 IRS source using the IRAF 
task BACKGROUND. 

In the H-band spectral segments, only the [Fe n] 1 .5999 and 
1.6440yum lines have been detected, while several additional 
transitions have been observed in K-band. Figures [T] and [2] re- 
port the K-band spectral images showing the different emission 
lines found for the HH34 and HH1 jets, while the most promi- 
nent identified lines are listed in Table [T] for both the jets. For 
HH34, we report here only the lines which are spatially resolved 
along the jet direction: additional transitions, mostly from per- 
mitted species, have been detected in the on-source spectrum and 
discussed in Antoniucci et al. ( |2008[ ). Along the HH34 jet, our 
detected lines are consistent with those observed by Takami et al. 
( 2006 ) except for our additional detection of the [Fe n] 2.224 [im 
line, originating from the 2 H term. The [Fen] 2. 133 /mi line, 
connecting the 2 P and 4 P terms, has been detected in both the 



HH34 and HH1 jets: as pointed out by Takami et al.|([2006), this 



transition, with excitation energies in excess of 25 000 K, is very 
useful to probe high excitation regions, such as those found at 
the jet base or in high velocity shock interaction. In HH1, we 



1 IRAF (Image Reduction and Analysis Facility) is distributed by 
the National Optical Astronomy Observatories, which are operated 
by AURA, Inc., cooperative agreement with the National Science 
Foundation. 



also detected two transitions from [Ti n], characterised by exci- 
tation temperatures on the order of 7000 K. Finally, in knot H 
of HH1, weak emission from Bry is detected, testifying for the 
high excitation conditions in this knot. 

3. [Fen] and H 2 kinematics 

In order to study the kinematics of the two jets in both 
the atomic and molecular components, we have constructed 
position- velocity (PV) diagrams of the [Fen] 1.644/im and H2 
2.122/im emission lines. The velocity is expressed with respect 
to the local standard of rest for both PV diagrams. A parental 
cloud velocity of 8km s _1 and!0.6kms _1 has been adopted for 



HH34 and HH1, respectively ( An glada et al.|1995||Choi & Zhou 
|1997[ ), and subtracted in the final PV velocity scale. Distance 
scales (in arcsec) have been measured with respect to HH34 
IRS and VLA1 for HH34 and HH1. Since the driving source 
VLA1 of the HH1 jet is not visible in the near-IR, we have used 
the bright Cohen-Schwartz (CS) star as the positional reference, 
adopting for it an angular separation of 35'.'9 from the source 
VLA1. In the following we will describe the results obtained for 
the two sources separately. 

3.1. HH34 



3.1.1. Large scale properties 

Figure [3] and |4] show the [Fen] 1.644 /on and H2 2.122yum lines 
PV diagrams relative to the blue-shifted and red-shifted lobes of 
HH34. The intensity scale is different in the two figures to evi- 
dence the relatively weak red-shifted emission. The redshifted 
counterpart of HH34 jet is clearly detected in our spectral im- 
ages. We have named the detected blue-shifted emission knots 
from A to L, following the nomenclature of Eisloffelj&Mundt] 
(1992J and Reipurth et al. (2002). We do not detect emission, 
however, from the knot B, since this knot is not aligned with the 
main jet axis, as evidenced in Reipur th" et al.| ([2002). The next 
knots we observe are the knots C and D that we have group to- 
gether as knot C. The red-shifted knots, which are approximately 
located at symmetric positions with respect to the corresponding 
blue-shifted knot, are named here from rA to rH. 

The values of the [Fen] and H2 1-0S(1) peak velocities have 
been computed applying a Gaussian fit to the line profile of ev- 
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Table 2. Observed radial velocities along the HH34 jet. 



Blue lobe Red lobe 







[Fe n] 1 .64 /jm 


H 2 2.12 








[Fen] 1.64 //m 


H 2 2.12//m 


Knot 


rFl 


HV H_J 


T XI 0" 


HVC* 


]_S/c h 


Knot 


a 
L t 




hvc* i vc b 

1 I V ^ . J i V V - 













-A 


rA 


+ 1 




-i-10 


A3 


-3.0 


-98 (37) 


-67 


-98 (17) 


-15 


rB 


+7.3 


+ 141 (51) 




Al 


-4.5 


-93 (37) 


-68 


-95 (12) 


-7 


rC 


+ 12.3 


+ 128 (58) 




C 


-9.0 


-92 (43) 




-94 (12) 


-7 


i-D 


+15.0 


+ 108 (35) 




E 


-12.5 


-108 (44) 




-110(8) 




i-E 


+ 17.4 


+96 (43) 




F 


-15.2 


-100 (43) 




-107 (12) 




rFG 


+21.0 


+96 (46) 


+ 113 (26) 


G 


-17.2 


-106 (43) 




-111(8) 




rH 


+25.4 


+ 100 (35) 


+ 115 (21) 


I 


-20.0 


-100 (43) 




-105 (24) 












J 


-22.1 


-96 (43) 




-99 (19) 












K 


-24.8 


-92 (44) 




-92 (26) 












L 


-28.9 


-98 (26) 




-99 (19) 













" Distance from the source in arcsec given by the mean value in the adopted aperture. Negative values correspond to the southeastern, blue- 
shifted jet axis. 

* Radial velocities (in km s _I ) with respect to the local standard of rest and corrected for a cloud velocity of 8 km s~ ! . The radial velocity error 
is 2kms~' . LVC and HVC refer to the Low Velocity Component and High Velocity Component. The velocity dispersion (in kms~' ), measured 
from the line FWHM deconvolved for the instrumental profile, is reported in brackets for the HVC. 



ery blue and red-shifted knot emission, and they are listed in 
Table [2] In the knots closer to the star, where two velocity com- 
ponents have been identified, their peak velocity has been sepa- 
rately measured considering a two-Gaussian fit. The [Fe n] radial 
velocities in the blue lobe cover a range from -92 to -108 km s -1 
, which is consistent with the values measured by Takami et al. 
(2006| and Davis et aL] ( |2003| l, who found a range in radial ve- 
locities for the blue lobe from -90 to -lOOkms -1 , corrected for 
a cloud velocity of 8 km s . 

From knots A6 to I the blueshifted radial velocity increases 
from — 92 to — lOOkms -1 , passing through a maximum at 
~ -108kms~' , then decreases again down to — 92kms~ 1 at 
knot K. Finally it increases in knot L to — 98 km s _1 . Errors in 
relative velocities are estimated to be on the order of 2kms _1 , 
since our wavelength calibration has an uncertainty of 0.1 A. 
The [Fen] radial velocity along the red lobe shows a similar 
behaviour. The velocity roughly increases from the knot clo- 
sest to the source, rA, to the knot rB, from ~130kms _1 to 
~140kms _1 , then decreases down to a value of ~96kms~' . 
Table 2 reports the velocity dispersion of the [Fe n] emission, 
measured from the FWHM of the Gaussian fit, deconvolved 
for the instrumental profile. In knots A6-A3-A1, where diffe- 
rent velocity components are evident (see Section 3.1.2) the re- 
ported velocity dispersion refers to the brightest component at 
high velocity. Intrinsic line widths on the order of 35-40kms _1 
are observed all along the jet. Under the assumption that the 
line emission arises from unresolved shock working surfaces, 
we have, following Hartigan et al] ( |1987| l, that the shock velo- 
city is roughly given by V s ~ AV(FWZI)~ 2 x AV(FWHM). 
This implies shock velocities of the order of 70-80km s~ 1 , thus 
much higher than the value of ~30km s _1 estimated by |Hartigan| 
et al. ( 19941 on the basis of the comparison of optical line ra- 



tios with shock models. Indeed, shocks with speeds as high as 
80km s are expected to produce a strong ionisation, on the or- 
derof x e ~ 0.3 - 0.4 (Harti gan et al.||1994| ) while |Podio et al. 
2006 ) measured an average ionisation along the HH34 jet of 



mined not only by the shock but also by, e.g. a lateral expansion 
of the jet. 

H 2 radial velocities have a range from -89 to - 1 1 km s~ 1 for 
the blue lobe and an average value of — h 1 15 km s _1 for the red 
lobe. The H 2 line profile is resolved all along the jet, with AV 
values of the order of 10-30 km s . 

Also in H 2 , radial velocities show cyclic variations on small 
scales along the jet, with an increase of roughly 20 km s from 
knot A6 to knot E and a subsequent decrease of the same order 
at knot K. It is known that the HH34 jet presents velocity vari- 



ability on large and small spatial scales. Raga & Noriega-Crespo 



(1998) and Raga et al. (2002) have shown that to reproduce the 
velocity pattern observed at different epochs, a model of vari- 
able ejection velocity including three modes with different pe- 
riods is needed. The fastest of these modes can be represented 
with a sinusoid having a period of 27 yrs and an amplitude of 
~15kms _1 . This is roughly consistent with our observed velo- 
city pattern, that is reproduced quite closely by both [Fe 11] and 
H 2 and by the red-shifted and blue-shifted gas, clearly indicating 
its origin from ejection velocity variability. The velocity pattern 
observed at large distance may, however, be biased by the change 
in the jet axis direction that occurs at d~5" from the source, cou- 
pled with the widening of the jet diameter (up to ~076 at ~20 



arcsec from the source, Reipurth et al. 

or 2 



strumental slit of only 073, part of the 



|2002 1. Adopting an in- 



cinematical components 



only 0.04. Therefore, it seems that the line widening is deter- 



of the jet may not be properly probed by our observations. 



3.1.2. Small scale properties 

Close to the central source, the [Fen] lines broaden and emis- 
sion at lower velocities, down to Okms -1 appears within ~3" 
from the central source (Figure BV Inside ~\" we see emis- 
sion also at positive velocities, reconnecting with the spatially 
resolved red-shifted knot rA. This central region was already ob- 
served in_[Fe n] and H 2 by |Davis et aT] ( [200T][2003] i and |Takami| 
et al. (2006l. Our observations have, however, a better spatial 
resolution than Davis et al. (2001 2003) and are much deeper 
than those of Takami et al. (2006 ), who did not detect the [Fe 11] 
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Fig. 3. Continuum-subtracted PV diagrams of the [Fen] 
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Fig. 4. Continuum-subtracted PV diagram for the H2 1-0S(1) 
emission line for the blue and red lobe of the HH34 jet. A P.A. 
of -15° was adopted. Contours show values of 5, 15, 45, 135, 
405 cr for both the lobes. On the Y-axis the distance from HH34 
IRS is reported. 



redshifted component that we see at the jet base. We identify 
the high blue-shifted velocity corresponding to the large scale 
jet as the High Velocity Component (HVC) and the emission 
component from to ~50kms~' the Low Velocity Component 
(LVC), in analogy with the HVC and LVC observed in the for- 
bidden emission line (FEL) regions of T Tauri stars. We remark, 
however, that in T Tauri FELs one usually observes two sepa- 
rated peaks spatially located at different offsets from the central 
source, the HVC being commonly displaced further downstream 
(e.g. Hirth et al. 1997 1. Here the HVC has a peak on-source while 
the LVC is seen as a weaker shoulder of the brightest component. 
This behaviour could, however, be due to the lower resolution 
of our observations as a consequence of the larger distance of 
HH34 with respect to other studied T Tauri stars. When com- 
pared with other [Fe n] PV diagrams observed in T Tauri stars, 



HH34 is more similar to HL Tau and RW Aurigae (Pyo et al. 
2006} than to L1551 and DG Tau ( |Pyo et al.|2002||2003) . 



At variance with [Fen], the H2 PV diagram shows spatially 
and kinematically separated LVC and HVC, and only the LVC 
is visible down to the central source. This component is close to 



0kms~' LSR velocity and the blue-shifted and red-shifted jets 
differ by less than lOkms 1 . The HVC appears at a distance of 
2" from the central source, at the position of the knots A6 and 
rA. At intermediate velocities between these two components, 
no emission is seen in the H2 PV diagram. This suggests that the 
two components correspond to physically distinct regions. The 
origin of the H2 high and low VC will be further discussed in 
Section 5.1. 



3.2. The HH1 jet 

Figures [6] and [7] show the PV diagrams of the H2 2.1 22 ymi 
and [Fen] 1.644 fxm lines of the HH1 jet. The PVs of the 
[Fe 11] 2.133 fj.m and [Tin] 2.160 ymi lines are also presented in 
Figure [7] The observed knots are named from KL to P follow- 
ing the nomenclature by Eisloffel et al. ( 1994| l and Bally et al. 
(2002). [Fen] is detected only in the knots closer to the star, 
from KL to E, while H2 can be traced all along the jet; as al- 
ready shown by Davis et al. ( 2000[ > and Nisini et al. ( 2005 1, the 
ratio [Fen] /H2 sharply decreases with the distance from VLA1. 
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the region nearest to the source. A P.A. of -15° was adopted. 
Contours show values of 5, 15, 45, 135, 260 <x for both lines. On 
the Y-axis the distance from HH34 IRS is reported. 



Table 3. Observed radial velocities along the HH1 jet. 



Knot 


[Fe n] 1 .64 fim 

vfl 


H2 2.12 fim 
Blue Red 


KL 


1.6 


-28 (64) 


-16 


+ 13 


JI 


3.4 


-46 (53) 


-27 


+6 


H 


5.1 


-43 (47) 


-24 


+ 14 


G 


7.5 


-35 (22) 


-17 


+23 


F 


9.8 


-31 (20) 


-14 


+ 17 


E 


11.6 


-24 (29) 


-11 




D 


13.6 




-12 




C 


15.0 




-10 




B 


17.0 




-11 




A 


18.2 




-3 




M 


22.1 




-4 




N 


23.5 




-12 





" Distance from the source in arcsec given by the mean value in the 
adopted aperture. 

6 Radial velocities (inkms -1 ) with respect to the LSR and corrected 
for a cloud velocity of 10.6 km s . The radial velocity error is 2 km s _I . 
The velocity dispersion (in kms~' ), measured from the line FWHM 
deconvolved for the instrumental profile, is reported in brackets for the 
HVC. 



Such a decrease is accompanied by a decrease of the H2 and 
[Fen] radial velocities. Also the velocity dispersion diminishes 
from the internal to the external knots, as shown in Table [3] In 
particular, we measure intrinsic FWHM (i.e. deconvolved by the 
instrumental width) decreasing from ~ 64 km s _1 in knots KL, to 
~ 20kms _1 in knots G-F. As in HH34, these values imply high 
shock speeds (from ~ 120kms _1 in knot KL to ~ 40kms~' in 
knot G), progressively decreasing outwards. Indications of shock 
velocities in the HH1 jet larger than in HH34 for the internal 
knots are given by the detection of the Bry and [Fe 11] 2.133 fim 
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Fig. 6. Continuum-subtracted PV diagram for the H2 1-0 S(l) 
line along the HH1 jet. A position angle of 145° was adopted. 
Contours show values of 10, 40, 60, 100 cr. On the Y-axis the 
distance from the driving source VLA1 is reported. 



lines in knot H. These lines are not detected in knots KL and 
JI probably due to the larger extinction. The [Fe 11] and H2 radial 
velocities appear double-peaked from knot KL to H in [Fe n] and 
from knot KL to F in H2 . Both, [Fe 11] and H2 present one of the 
two components red-shifted with respect to the LSR. 

The [Fe 11] blue-shifted radial velocity increases from knot 
KL to knot JI from a value of — 28kms~ 1 to — 46kms _1 and 
then decreases again to a velocity of — 24km s _! in knot E. The 
red-shifted component is, however, too faint to be fitted with a 
Gaussian profile. Nevertheless, we can estimate an average value 
of ~ +30 km s -1 from knot KL to H. 

On the other hand, in the H2 blue-shifted component the 
velocity increases from knot KL to knot JI, as in the case of 
[Fe 11] and then roughly decreases to — 3 km s _1 in knot A. After 
that, the radial velocity increases again to a value of around 
-12kms _1 in knot N. The red-shifted component covers a ve- 
locity range from +6 to +23kms~ 1 . The H2 line profiles are 
represented in Figure [8] for several knots in the region closest 
to the source. The measured H2 velocities are consistent with 
the ones reported in |Davis et"aL ( 2000 1 for the knots F and A. 
The H2 lines are non-resolved implying an intrinsic line width 
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Fig. 7. PV diagrams for the [Fen] 1.64 /mi, 2.13 fim and [Tin] 2.16/mi emission lines of the HH1 jet. A position angle of 145° 
was adopted. Contours show values of 11, 22,. .., 176 cr; 4, 8, 12cr and 3, 6, 12cr, respectively. On the Y-axis the distance from the 
driving source VLA1 is reported. 



less than 39 km s 1 . The origin of the double velocity compo- perpendicular to the line of sight. These two possibilities will be 

further discussed in Section 4.1.2. 
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Fig. 8. Normalised line profiles of the H2 1-0S(1) for different 
knots near the source. 



nent can be due either to the presence of another jet (e.g. the one 
driving the HH501 object |Reipurth et al.|2000{|Baliy et al.|2002) 
or to the emission along the wings of a bow shock seen almost 



4. Diagnostics of physical parameters 

4.1. Electron density 

The electron density in the atomic jet component can be derived 
from the ratio of the [Fen] 1.600/1.644 /mi lines. This ratio is 
sensitive to n e values between ~10 3 and 10 s cm" 3 , while it de- 
pends only weakly on the temperature (e.g. Nisini et al.|2002 1. 
Electron densities of the [Fen] emission line region as a func- 
tion of the distance from the driving source has been measured 
by |Podio et aL1 ( (2006] l and|Nisini et al.|p005) for the HH34 and 
HH1 jets. In addition, Takami et al.|(|2006) provide the electron 



density in the inner region of HH34. Taking advantage of the 
velocity resolved profiles in both the 1 .644 and 1 .600 fim lines, 
we can now measure the electron density in the different velo- 
city components. To do that, we have extracted the spectra of the 
two lines at different positions along the flows and measured the 
1.600 /mi/1.644 /mi intensity ratio in each pixel along the spec- 
tral profile. Figures|9]and[T0]show the normalised profiles of the 
two lines and their ratio as a function of velocity for HH34 and 
HH1. The spatial intervals used to extract the spectra of the in- 
dividual knots are given in Tables |4] and [5] The [Fen] line ratio 
has been computed only for the velocity points where the inten- 
sity in both the lines has been measured with a S/N larger than 
three. The plotted ratio gives a qualitative indication on how n e 
varies in the different velocity components, that is, a higher ra- 
tio indicates a higher electron density. In the internal knots of 
HH34 (from A6 to A3) the 1.600 /im/1.644 fim ratio decreases 
by ~70% going from — 50kms~' to -100 kms" 1 . In the knots 
at larger distance from the central source, on the other hand, the 
maximum 1.600 /mi/1.644 fim ratio is observed at the radial ve- 
locity peak, with some evidence that the ratio decreases in the 
line wings at both higher and lower velocities. In HH1, the 1.600 
fim line has been detected with a S/N larger than three only in the 
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Fig. 9. Normalised line profiles (lower panels) of the [Fe n] lines 1 .644 fim (solid line) and 1 .600 /im (dotted line), and their ratio in 
each velocity channel (upper panel) for different extracted knots along the HH34 jet. Note that rA is a region that includes the knot 
rA. 



internal knots, from JI to F. Here, the 1 .600 /im/1 .644 //m ratio 
has a minimum close to Okms -1 velocity and increases, up to a 
factor of two in knot KL, towards high velocities. This behaviour 
is therefore different from what has been observed in the inter- 
nal knots of HH34. In HH1, however, the inner jet region is not 
detected due to the high extinction and therefore we are not ob- 
serving here regions at the jet base as in HH34. It seems indeed 
that the dependence of the electron density with the velocity is 
different between the FEL regions close to the star and the knots 
along the jet beam. 

In order to have a quantitative determination of the electron 
density in the different velocity components, we have separated 
in the line profile the contribution from the HV and the LV com- 
ponents, and applied a 16 level Fe + statistical equilibrium model 
( Nisini et al.||2002] l. In the following, the results are discussed 
separately for the two sources. 



— 7kms -1 (LV). The adopted A v and T e values are given in 
Table |4j together with the derived electron densities in the two 
components. Figure [TT] plots the derived values of n e as a func- 
tion of the distance from the source. The first trend that we no- 
tice is a sharp decrease of the HVC electron density (from ~10 4 
to ~ 2 x 10 3 cm -3 ) from the knot A6 to the other knots at dis- 
tances farther than 2'.' 5 from the source. Such a decrease in n e has 
also been observed in Podio et al. ( 2006[ > and the values they de- 
rived, scaled for the different considered spatial regions observed 
through slits of slightly different width, are consistent with our 
measured values for n e . In the red-shifted knot rA, which is the 
only knot where significant 1.600 /urn emission has been detected 
in the red-lobe, we find a value for n e ~ 1 .2 x 10 4 cirT 3 , i.e. com- 
parable to the value derived in knot A6. Secondly, as shown in 
Figure 



HVC: 



11 



the values for n e in the LVC are higher than in the 
in knot A6, the LVC electron density is 2.2 x 10 4 cm" 3 , 
i.e. 70% higher than in the HVC. About the same percentage is 
measured in knots Al and A3. From the information in our data, 
we are unable to disentangle, whether the larger electron den- 
sity in the LVC with respect to the HVC is due to a higher total 
density or to a higher ionisation fraction. Our result is, however, 
in agreement with what is found in the FEL regions of T Tauri 
stars, where the emission component associated with the low ve- 
locity gas is denser and less excited than the HVC (Hamann et al. 
1994; Hartiga n et al.|1995) . In such studies, the LVC was not 
spatially resolved and the different values for n e found between 
-120km s _1 to — 66kms~ 1 (HV) and from — 66kms~' to the HVC and LVC were interpreted assuming that the LVC orig- 



4.1.1. The HH34jet 

In HH34, we have defined the HVC and LVC velocity ranges 
from the profile of the 1.600 jjm line of knot A6, where the 
two components have been fitted with a two-Gaussian fit. We 
have then measured the line ratios considering the intensities 
integrated in the FWHM ranges of these two Gaussians also 
for all the other knots. The considered velocity bins are from 
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Fig. 10. Normalised line profiles (lower panels) of the [Fe 11] lines 1 .644 /urn (solid line) and 1 .600 /mi (dotted line), and their ratio 
in each velocity channel (upper panel) for different extracted knots along the HH1 jet. 



inates from a dense compact region close to the disc surface, 
while the HVC is associated with a more extended high velocity 
jet displaced further out. Such an interpretation was supported by 
the spatial offset from the central source often found between the 
two components of the forbidden lines detected in spatially unre- 
solved spectra of T Tauri stars. Spatially resolved measurements 
of the electron density in the HV and LV components have been 
performed in DG Tau by Bacciotti et al. ([2000) and |Lavalley- 



|Fouquet et al.| l |2000) . In this case, at variance with HH34, the 
electron density has been found to increase with velocity, up to 
a distance of ~ 450 AU from the central source. On the other 
hand, in the DG Tau micro-jet the total density in the LV com- 
ponent is higher than in the HV component, because of the much 
lower ionisation fraction. Therefore, it might be that the density 
structure in HH34 and in the T Tauri stars is similar, but they 
differ in the excitation conditions. In HH34, we spatially resolve 
the LVC, whose emission is traceable up to ~5 arcsec from the 
source, i.e. more than 2000 AU. Therefore, the spatial scale be- 
tween the DG Tau micro-jet and the HH34 jet are very different: 
in addition, and at variance with many T Tauri stars, we do not 
observe a spatial offset between the HV and LV components, al- 
though this may result from the moderate spatial resolution of 
our observations. Therefore, it may be that the LVC that we are 
tracing with our observations has a different origin than in the 
CTTS case. In the outer knots, i.e. from C outwards, we can- 
not distinguish High and Low velocity components anymore. 
Instead, the density here seems to have an opposite behaviour, 



with the higher density at the velocity peak and the lower den- 
sities in the line wings. The analysis performed in these knots 
located far from the source could, however, be affected by a not- 
perfect alignment of the slit with the jet axis and by the intrin- 
sic jet width larger than the slit, as discussed in Sect. 3.1.1 . 
Nevertheless, such a pattern agrees with the results obtained by 
Beck et al. ( 20Q7\ using integral field spectroscopy of this part 
of the jet. They found that both the velocity and the electron den- 
sity decrease with distance from the jet axis. Such a behaviour is 
consistent with models for jet internal working surfaces. 



4.1.2. The HH1 jet 

In HH1 we define a low and high velocity component from the 
[Fe 11] 1 .644 /im line profile of knot KL. We selected a HVC 
from — 80 to — 22 km s _1 and a LVC from ~-22 to ~36 km s _1 . 
In analogy to HH34, n e for the LVC and the HVC has been com- 
puted separately for all extracted knots. Visual extinction and 
electron temperature values in each knot have been taken from 
Nisi ni et ak] ( |2005[ ). At variance with the HH34 jet, the n e in the 
LVC is lower than in the HVC. In the HVC n e decreases knot by 
knot from a value of 10 4 cirT 3 in knot KL to 3.8 x 10 3 crrT 3 
in knot F. On the other hand, in the LVC n e decreases from 
9.8 x 10 3 cm" 3 to 8.3 x 10 3 cm 3 in knot KL and JI, respectively. 
These n e -values are in agreement with those derived by |Nisini 
et al. (2005) in velocity integrated spectra. 
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Table 4. M jet along the HH34 jet. 



Knot 


r,Fl 


A,F1 


V.fHVcFl 

(kms -1 ) 


n fHVCFl 


m .nwcfl 

lvl jelV 11 v -M 

(Mo vr~') 


v.itvct 


n (iyo d 








map 


(10 3 cm"') 


(kms -1 ) 


(10 3 crrT 3 ) 


(M„ vr^'l 


A6 


(-2.6, +0.9) 


7.1 


218 


10.5 


5.2 x 10~ 8 


124 


22.5 


6.6 x 10" 9 


A3 


(-3.8,-2.6) 


7.1 


232 


1.8 


8.3 x 10~ 8 


159 


5.0 


2.0 x lO -9 


Al 


(-5.4,-3.8) 


7.1 


220 


1.8 


3.3 x 10~ 8 


162 


5.0 


1.3 x 10~ 9 


C 


(-9.9,-8.0) 


7.1 


218 


1.2 


4.4 x 10~ 8 








E 


(-14.1,-11.0) 


1.3 


255 


1.8 


5.9 x 10~ 8 








F 


(-16.6,-14.1) 


1.3 


236 


1.8 


5.0 x 10~ 8 








G 


(-17.8,-16.6) 


1.3 


251 


2.5 


1.0 x 10~ 7 








I 


(-21.0,-19.2) 


1.3 


236 


2.5 


4.2 x 10~ 8 








J 


(-23.3,-21.0) 


1.3 


227 


1.2 


3.9 x 10~ 8 








K 


(-26.1,-23.5) 


1.3 


217 


1.5 


l.Ox 10~ 8 








L 


(-30.3,-27.5) 




232 












rA 


(0, +2.0) 


7.1 


218 


12.0 


2.2 x 10~ 9 









Distance from the sou rce in arcsec. Neg ative values correspond to the southeastern, blue-shifted jet axis. 



Visual extinction from 



Podio et al 



(2006) 



Tangential velocity assuming an inclination of the jet i=22°7 to the plane of the sky (Eisloffel & Mundt 
d Electron density for the HV and LV components. 

' Mje, for the HV and LV components assuming an electron temperature of 7000 K. 



1992). 



The presence of two velocity components and the corres- 
ponding dependence on electron density with the velocity, was 
already found by |Solf et aL| ( 1991) , who produced a PV dia- 
gram of the electron density in the HH1 jet using the optical 
[S n] A 6716/6731 line ratio. They found a blue-shifted compo- 
nent with an electron density of ~4000 crrT 3 and a slightly red- 
shifted component having a lower density of ~ 1000-2000 cm -3 . 
They interpreted the blue-shifted high-density component as due 
to scattered light originating from a jet region closer to the star, 
where the density and excitation are higher. This interpretation 
is, however, difficult to keep in the light of our IR observations. 
Indeed, the blue-shifted line component is clearly identifiable 
as the main jet component, extending at large distances from 
the exciting source. Moreover, our H2 PV diagram (Figure [6} 
clearly shows the presence of two separate components that can- 
not be attributed to scattered light contribution. The redshifted 
component peaks at the KL position and decreases in intensity 
further out, while the second blue-shifted component gradually 
increases its intensity with the distance. One possible interpreta- 
tion is that the red-shifted component belongs to a different jet 
that intersects the HH1 jet at the KL position. Such a jet could 
be responsible for the two bright knots designated as HH501 ob- 
jects by Reipu rfhet al.|p000) , and that are moving away from 
VLA1 with a proper motion vector inclined with respect to that 
of the HH1 jet by ~10°. If we take the axis of this separate jet 
equal to the direction of the proper motion vector, the HH501 jet 
should intersect the HH1 jet at the position of the KL knots. 



4.2. Mass flux 

The mass flux rate, Mj et , along the beam of the HH34 and HH1 
jets has been recently measured in Pod io et aL]p006) andlNisini 
et al. ( 2005 1, respectively, using different tracers, both optical 
([S11] , [Oi] ) and infrared ([Fen] , H2). These works show that 



the mass flux derived directly from the [Fen] line luminosity, 
using the measured tangential velocity, is always equal or larger 
than the Mj el value derived from the luminosity of the optical 
atomic tracers, in spite of the possibility that part of the iron is 



still locked on dust grains. This is probably due to the fact that 
[Fe 11] traces a larger fraction of the total flowing mass than the 
optical lines, as discussed in Nisini et al. (2005 ). At the same 
time, it was found that in these jets the mass flux traced by the H2 
molecular component is negligible with respect to the mass flux 
due to the atomic component. Taking advantage of the velocity 
resolved observations available, we now want to measure Mj el 
in the different velocity components and examine which of them 
is transporting more mass in the jet. 

Mj et have been obtained here from the luminosity of the 
[Fen] 1.644 /im line, adopting the relationship M = firriH x 
(«hV) x v,/l t , where fi = 1.24 is the average atomic weight, 
niH and n# are the proton mass and the total density, V is the 
emitting region and v, and I, are the velocity and length of the 
knot, projected perpendicular to the line of sight. The product 

n H V can be expressed as L(line) (hvAifi ^ ^) . A,-, f, are the 
radiative rates and fractional population of the upper level of the 
considered transition and is the ionisation fraction of the iron 

Fe 

having a total abundance with respect to hydrogen of ^J. 

We have assumed th at all iron is ionised , and has a solar 
abundance of 2.8 x 10" 5 ( |Asplund etal.|2005| i.e. no dust deple- 
tion). This latter hypothesis leads to a lower limit of the actual 
mass flux, since it has been shown that the -velocity averaged- 
gas-phase iron abundance in the HH34 and HH1 jets might be 
only between 20 and 70% of the solar value ( Nisini et al.|2 005 ; 
Podi o et al.|2006) . We do not adopt the Fe abundance estimate 
provided in these papers, since we do not know if the depletion 
pattern is constant among the low and high velocity components. 

To compute the fractional population, we have taken the n e 
values derived separately for the LV and HV components, while 
we have assumed a single temperature for both components, 
equal to the values derived in |Nisini et al. (2005 1 for the HH1 
jet and a constant value of 7000 K for HH34 given by the av- 
erage value derived by Pod io et ah] ( [2006 ) in knots from E to 
I. 

Tangential velocity values have been derived from the ra- 
dial velocities assuming an inclination angle i=10° for HH1 



Table 5. M jet along the HH1 jet. 
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Knot 


rFl 

<H 


U 

mag 


T b 
(10 3 K) 


V f (HVCf] 

(kms-'F 


n e (HVCfl 

c V " A 1 

(10 3 cm"') 


M, e ,(HVCfl 

j el x A l 

( M yr 1 ) 


V,(LVC) e 
(kms- 1 ) 


n e (LVC)) d 
(10 3 cm" 3 ) 


MjJLVCY 
(Mo yr 1 ) 


KL 


(0.7,2.6) 


8.3 


11 


159 


10.6 


1.3xl0 -8 


187 


9.8 


2.3xl0- 9 


JI 


(2.6,4.3) 


8.3 


11 


261 


9.0 


3.6xl0~ 8 


187 


8.3 


3.1xl0" 9 


H 


(4.2,6.1) 


2.9 


9.2 


244 


8.3 


4.4xl0~ 8 








G 


(6.1,8.9) 


2.0 


10.5 


198 


4.6 


2.7xl0~ 8 








F 


(8.9,10.7) 


2.0 


9.8 


176 


3.8 


2.2xl0~ 9 









Nisini et al. 



(2005) 



" Distance from VLA1 in arcsec. 
* Visual extinction and electron temperature from 
c Tangential velocity assuming an inclination of the jet 1=10° to the plane of the sky (Bally et al. 
d Electron density for the HV and LV components. 
' Mjc for the HV and LV components. 



2002) 



and i= 22?7+5° for HH34 ( |Bally et aL|2002||Eisloffel&M undt 
1992|l. In the case of HH34, the inclination angle derived by 



Eisloffel & Mundt|([l9"9"2"l) is prefer red to the value of i=30° es- 



timated by Heathcote & Reipurth ( 1992), who did not consider 
the pattern motion of the jet knots. The resulting tangential ve- 
locities of the HVC derived here are listed in Table 4. Compared 
to the values derived by the proper motion analysis of the [S n] 
emission by Reipurth et al. ( 2002 ), we derive smaller velocities 
in the internal knots (by ~ 20-50kms 1 in knots A6-A3) and 
larger in the more distant knots. In particular, we see no evi- 
dence for a deceleration of the jet as that shown by Reipurth 



et al. (2002 1, and the velocity remains higher than 200km s 



along the jet length. This could be an indication either that the 
inclination angle does not remain constant along the jet, or that 
the [Fe n] line has a kinematical behaviour different from [S n] . 

The luminosity of the line has been computed by integrating 
the extinction corrected flux of the knot in the same range of 
velocity used to calculate the electronic density. In Tables |4] and 
[5] the derived values, together with the parameters adopted, are 
reported for the HH34 and HH1 jets. In Figure [TT] the HH34 
Mjetis plotted for the High and Low velocity components as a 
function of the distance from the central source. 

In both jets, M Jet (LVC) is lower than Mj et (HVC) by a factor 
of about 6 (in HH1) and 10 (inHH34). In HH34, the M jet (HVC)/ 
Mj e t (LVC) ratio does not reflect the derived n e ratio in the two 
components, indicating that the (n//V) x v, product is signifi- 
cantly smaller in the LVC (around 8 times) than in the HVC. 
Since the difference in tangential velocity between each compo- 
nent is not enough to justify such a result, the larger Mj et (HVC) 
indicates a significant higher {n H V) product in the HVC, i.e, the 
HVC has either an emission volume or a total density larger than 
the LVC, or both. This result could be biased by the finite slit 
width with respect to the jet diameter. In fact, the jet width mea- 
sured only by HST in the optical lines is less than 073 up to a 
distance of 5"from the source ( |Reipurth et al.||2002) . In addi- 
tion, the jet emission in the low velocity component could be 
broader than the jet width derived from velocity integrated emis- 
sion maps as observed, e.g. in DG Tau ( Bacciotti et al.| 2000). 
This effect propably does not account for the difference of nearly 
two orders of magnitude in derived mass-loss rates fro the LVC 
and HVC in the inner regions of HH34 but the derived mass-loss 
rates for the low-velocity component are likely lower limits. We 
also note that the reported values of My er (HVC) remain roughly 
constant (~ 5 x 10~ 8 M /yr) along the whole jet as expected in 
steady jet flows. This fact excludes significant flux losses as the 
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Fig. 11. The electron density (upper panel) and mass flux (bot- 
tom panel) are represented as a function of the distance from the 
source for the HH34 jet. Solid circles indicate the the electron 
density and mass flux for the HVC, while open circles refer to 
the values of the LVC. 



jet opening angle increases. In jets from CTTSs, Hartigan et al. 
( 1995 ) have already shown that the mass flux in the LVC should 
be lower than in the HVC: from the analysis of not spatially re- 
solved optical spectra, they concluded that the LVC in jets from 
CTT stars could be responsible for carrying the majority of mass 
and momentum only if the emitting region were smaller than 1 
AU. 

In HH34, the mass flux derived for the internal knot A6 is 
very low, much smaller than the Mj et value in knot A3. This is 
likely due to the fact that we assumed here the same extinction 
value of 7.1 mag that was estimated by Podio et al. ( |2006| ) as an 
average over the entire knot A, i.e. over 4". The extinction value 



on-source has been estimated to be about 45 mag ( jAntoniucci 
et al.||2008) : therefore a large dust column density gradient is 
expected in the inner jet region. Assuming an A v lower limit of 
7.1 mag also in the red-shifted component rA, we derive here a 
lower limit for Mj et of ~ 1 x 10~ 9 M yr -1 . The same problem 
probably applies to the internal knot KL in HH1, where the 8.3 
mag of extinction are the average value in the whole jet region 
LI. 
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We finally note that while the Mj et determined here for the 
HH34 jet agrees fairly well with the velocity integrated values 
estimated in |Podio et al. ( 2006| l, there is some disagreement be- 



tween the Mj et derived for HH1 here and in Nisini et al. (2005). 
The l arger discrepancy i s being found in the internal JI knots, 
where 
of the 



Nisini et al. 



( 2005 1 report a Mj et a factor of six larger. Part 
disagreement is due to the different adopted tangential ve- 
locity. The largest discrepancy is, however, due to a smaller flux 
(a factor of three) measured in the ISAAC Of.' 3 slit with respect 
to the 1" slit used in |Nisini et al.| ( |2005] l. It can be that this latter 
measurement was contaminated by the presence of the second 
jet responsible for the redshifted velocity component or that the 
HH1 jet itself has a diameter wider than 073. 

It is interesting to compare the derived mass flux rates with 
the source mass accretion rate, in order to estim ate the M / e , /M acc 
efficiency in embedded young sources. Antoniucci et al. (2008), 
have derived a mass accretion rate of the order of 4.1 x 10~ 6 
M yr 1 for HH34 IRS. This would imply M jet /M acc > 0.01, 
in agreement with what is found in T Tauri stars and predicted 
by MHD jet launching models (e.g. [Ferreira et al.|2006 1. There 
are no measurements of the mass accretion rate in VLA1, the 
exciting source of the HH1 jet, available. We can give a rough 
estimate of this value by assuming that the source bolomet- 
ric luminosity (L/, / ~ 50 Lq, |Chini et~aL 1997\ is entirely 
due to accretion. Such an assumption is supported by the fact 
that VLA1 is a known class source. If we further assume a 
stellar mass and radius of 1 Mq and 3 R© respectively, we 
get M acc ~ 6 10~ 6 M Q yr~', which would imply in this case 
Mjet I > 0.007. 



4.3. The Ti + /Fe + ratio in HH1 

Forbidden [Ti n] emission lines have been detected for the first 
time in HH1. Similarly to iron, Ti has a low ionisation poten- 
tial, of only 6.82 eV, and thus is expected to be fully ionised 
in the jet plasma. The two detected lines have excitation ener- 
gies on the order of 7000 K, and critical densities, at 10 000 K, 
~ 5 x 10 4 cirT 3 . They have similar excitations, thus, with re- 
spect to [Fen] IR lines. Titanium, as iron, is a highly refractory 
element, thus the ratio of [Ti n] /[Fe n] lines can give some clue 
on the relative gas-phase abundances of Ti and Fe, and, in turn 
on the selective depletion of these elements on dust grains still 
present in the jet beam. Nisi ni et aE] ( 2005 ) have indeed shown 
that the gas-phase abundance of Fe in HH1 is lower than the 
solar value, in particular in the inner jet regions where it is only 
between 20%-30% of the solar Fe abundance. This indicates that 
part of the iron is depleted on grains and consequently that a sig- 
nificant fraction of dust is present in the jet beam. We can here 
check if titanium follows the same depletion pattern, by com- 
paring the observed [Tin] 2.160 //m/ [Fen] 1.644 fim ratio with 
the value theoretically expected assuming solar abundance va- 
lues given by |Asplund eTaT1 ( |2003] l, i.e. [Fe/Ti] Q = 354. For this 
analysis we have performed a statistical equilibrium calculation 
taking the [Ti n] radiative transition rates and electron collisional 
rates calculated in Bautista et al. (2006, Bautista, private com- 
munication), and considering the temperature and density values 
measured in each HH1 knot in Nisin i et al.| ( |2005| ). The [Tin] 
2.160 /im/[Feii] 1.644 /mi intensity ratio is rather insensitive to 
the adopted physical conditions, being on the order of ~ 500 in 
a temperature range from 9000 to 12 000 K, and density range 
from5xl0 3 tol.2xl0 4 cm . The observed ratios range instead 
between 150 (knot JI) and 280 (knot G), implying a gas-phase 
Fe/Ti ratio 2-3 times lower than solar. Thus there is an overabun- 
dance of Ti in the gas phase with respect to Fe relative to the solar 



value. This result indicates that the release of gas-phase elements 
from dust grains likely follows a selective pattern in which Ti- 
bearing condensates are more easily destroyed. A large Ti/Fe 
abundance ratio relative to solar abundances has been observed 
previously in the ejecta of rj Carinae by Bautista et al. (2006) 
who suggest two different scenarios to explain this finding: ei- 
ther there is a spatial separation between Ti- and Fe-bearing con- 
densates in the same dust grain, that makes the titanium more 
exposed to evaporation, or Ti-bearing grains are smaller than 
the Fe-bearing grains and thus more easily destroyed. Studies of 
this kind, employing a larger number of refractory species, seem 
therefore promising to constrain the composition and structure 
of dust grains in different environments. 

5. HH34 small scale jet: comparison with models 

The [Fe n] PV observed in HH34 is similar to the ones observed 



in some CTTSs, such as HL Tau and RW Aurigae (Pyo et al. 
2006). In these latter, as in other CTTSs observed in the op- 
tical, the HVC peak has an offset with respect to the central 
source, of typically 50-100 AU. Given the larger distance of 
HH34 with respect to the sources in Taurus, we are not able to 
resolve such spatial regions, and the HV peak in our PV is cen- 
tred on source. In CTT jets, the LVC is usually confined within 
d ^ 200 AU from the source. In HH34, on the contrary, the LVC 
at -60kms _1 persists at larger distances, up to 1000-2000 AU. 
In the following, we discuss these findings in the light of current 
theoretical models for the production of YSO jets. 

Magneto-centrifugal jet launching models, such as the disc- 
wind and X-wind models 



jFerreira||1997| |Shang et al.][T998 
Garcia et al.||2001 Pudritz et al.||2007| ) predict the presence of 
a broad velocity range close to the source position, correspon- 
ding to the uncollimated outer streamlines. In both these mod- 
els, however, the LVC is confined to a relatively small region at 
d % 200 AU from the source. Thus, the LV gas we observe at 
larger distance needs to be locally heated, probably by shocks. 

Within 400 AU from the HH34 IRS source the line emis- 
sion decreases with velocity almost continuously from its peak 

at lOOkms -1 to redshifted velocities up to — 1-200 km s -1 . 

We can compare the kinematical signatures we observe in 
this region with those in the synthetic PV diagrams available for 
the different jet models. Synthetic PV diagrams from X-wind 
models have been constructed for optical [Sn] and [Oi] lines 
( Shang et"aL"|l998|), wh ile [Fe n] PV plots have been presented 
by|Pesenti et al. \ 2003 i for the cold disc- wind model of Ferreira 
( |1997) . 

The cold disc-wind synthetic P-V diagram presented in 
Pesenti et al. (|2003[) was constructed assuming a spatial sam- 



pling and spectral resolution similar to those of our observations. 
This diagram predicts [Fen] 1.64 /an maximum (de-projected) 
velocities on the order of — 700kms _1 . In contrast, the maxi- 
mum velocity measured in HH34 from the FWZI in knot A6 is 
400 km s . However, the observed range of line of sight ve- 
locities predicted by disc-wind models significantly depends on 
the assumed range of launching radii and magnetic lever arm pa- 
rameter (A = (ta/Vo) 2 , with r A and r$ the the Alfven and launch 
radii, respectively). Ferreira et al. (2006J have shown that cold 
disc-wind models in general predict too large line of sight veloc- 
ities with respect to the observations, because they have too large 
magnetic lever arms of the order of 50 or higher. On the other 
hand, warm disc wind solutions with magnetic lever arms in the 
range 2-25 predict lower poloidal velocities, in better agreement 
with observations (Ferreira et al. 2006). In particular, our mea- 
sured poloidal velocities of the order of 250km s _1 , estimated 
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from the observed peak radial velocity corrected for the inclina- 
tion angle, would imply a value of A in the range 4-8, for laun- 
ching radii in the range 0.07-0.15 AU ( Ferreira et al.|2006 1. 



X-wind models predict a narrower spread in velocity that 
agrees with our observed range, taking into account the incli- 
nation angle of 23° of the HH34 flow. In such models, however, 
the intermediate and low velocity components are predicted to 
be present only in the inner ~ 30 AU from the central source, 
in contrast with observations. In addition, the synthetic PV di- 
agrams constructed by Shang et al. ( 1998] > do not include heat- 
ing processes in the calculation, but assume a uniform ionisation 
fraction and electron temperature. 

Another interesting feature predicted by X-wind models is 
the presence of a significant redshifted emission at the stellar 
position similar to what we observe in HH34. This is the first 
case where such a feature is observed, and indeed the absence of 
the redshifted emission in the PV of other CTT jets was taken 
as an indication to rule out the X-wind model for these sources 
(Taka mi et al.|2006| ). In HH34, such redshifted emission was not 
detected in previous observations probably because of the lack 
of sensitivity. 

In conclusion, a qualitative comparison based on the availa- 
ble synthetic PV diagrams indicates that both X-wind and disc- 
wind models reproduce the kinematics observed in the HH34 
jet at its base, while none of them can reproduce the large scale 
LVC observed at very large distances from the source. A more 
quantitative comparison would require the construction of [Fe n] 
synthetic P-V diagrams exploring a larger range of parameters 
and heating mechanisms. 

The other observational finding in our study that can be com- 
pared to model predictions, is the dependence of the electron 
density and mass flux on velocity. The fact that the electron den- 
sity is higher in the LVC than in the HVC is neither consistent 
with the available X-wind nor with disc-wind models. In disc- 
wind models the high velocity and dense gas ejected in the inner 
streamlines is bracketed by the slower gas at lower density orig- 
inating from regions of the disc more external and less excited. 
Cabrit et al. ( 1999) show, in fact, that the electron density in- 



creases toward the axis of the flow. Similarly, in X-wind models 
the electron density and ionisation decreases with the distance 
from the central axis ( |Shu et al7||1995| >. We point out that the 
main parameter predicted by different jet launching models is 
not the electron density, but the total density and its variation 
with velocity. Therefore, the found inconsistency may stem from 
an incorrect modelling of the excitation mechanism responsible 
for the gas ionisation, not necessarily from a wrong underlying 
MHD solution. 

An alternative possibility is that the LVC represents gas not 
directly ejected in the jet, but dense ambient gas entrained by 
the high velocity collimated jet. Observations at higher spatial 
resolution would be needed to perform an analysis similar to that 
presented here for the gas within 200 AU from the source. This 
would allow us to isolate as much as possible the component of 
gas ejected at low velocity in the jet itself. 



5.1. H 2 high and low velocity components in HH34 

High and low velocity H2 emission within a few hundred AU 
from the driving source has been detected in a dozen out- 
flows ( |Davis et al.||200T| ); such emission regions have been 
called Molecular Hydrogen Emission-line region (MHEL) in 
analogy to the atomic FELs regions observed in T Tauri stars. 
With respect to FEL regions, the gas components associated 



with MHELs usually show lower radial velocities and velocity 
spreads. 

In light of these previous findings, there are two scenarios 
that can explain the characteristics of the LVC in H2 lines in 
HH34. The first hypothesis is that such a low velocity warm mo- 
lecular gas is excited by oblique shocks occurring at the wall 
of a cavity created by the interaction of a wide angle wind with 
the ambient medium. Indeed, the large reflection nebulosity as- 
sociated with HH34 IRS suggests the presence of a cavity il- 
luminated by the central source. The second possibility is that 
the H2 gas originates from the external layers of a disc-wind. 
Magneto-centrifugal wind models in fact predict that the fast 
collimated inner component coexists with a slow and wide ex- 
ternal component at low excitation. In a previous study, the H2 
emission observed within 100 AU from DG Tau has been inter- 



preted by |Takami et al. (2004) in this scenario. Modelling of the 
H2 molecule survival and excitation in disc-winds is, however, 
needed to support this interpretation. 

The fact that the H2 HVC is observed only from a distance of 
~2" from the star, at variance with [Fe 11] , may indicate that the 
molecular gas does not survive in the conditions of the highly 
excited gas constituting the high velocity inner jet region. In this 
inner part the jet is probably travelling inside a low density cavity 
and thus no molecular material is present to be entrained by the 
fast moving flow. Reformation of H2 in the jet could be at the 
origin of the H2 HVC seen further out: we point out that dust 
in the jet is not fully destroyed (Podio et al. 2006) and thus the 
H2 molecules can be formed again once the physical conditions 
become favourable for the molecule survival. Excitation in the 
wings of mini -bow shocks along the flow can be another possible 
explanation for the H2 HVC, although the high radial velocity of 
H2 close to the velocity of the atomic gas, is not supported by 
H2 excitation models in bow-shocks ( |Flower et al.|2003"] ). 

6. Conclusions 

We have presented H and K-band spectra of the HH34 and 
HH1 jets, where [Fen] , H2 and for the first time [Tin] emis- 
sion lines have been detected. These observations provide us 
with detailed information about the kinematics of the emitting 
gas allowing us to resolve two velocity components in each jet 
and to measure their physical parameters separately. From the 
[Fe 11] 1 .600 fim /l .644 fim ratio the electron density along both 
jets and for each velocity component has been derived. In addi- 
tion, the mass flux has been inferred from the [Fen] 1.644 /mi 
line at different positions along the jet axis and for both velocity 
components. Finally, the [Tin] /[Fen] ratio gives us important 
information about the mechanism of dust reprocessing. Our re- 
sults can be summarised as follows: 

- In HH34, the atomic and molecular gas show two velocity 
components, the so-called high (HVC) and low (LVC) velo- 
city components, near to the source. We also detect for the 
first time the fainter red-shifted counterflow down to the cen- 
tral source and up to a distance of 25" from the star. 

- HH1 is traced down to ~ 1" from the source. The kinema- 
tics of the region closest to the driving source is again cha- 
racterised by a double velocity component, one blue-shifted 
and one red-shifted with respect to the source LSR. We in- 
terpret the red-shifted component as part of another jet, pre- 
viously detected by Reip urth et al.| ( |2000] ) in optical and near- 
IR lines. 

- In the innermost HH34 jet region, the electron density in- 
creases as the velocity of the jet decreases, with average 
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values of 10.8 x 10 3 cm- 3 and 4.7 x 10 3 cm- 3 for the LV 
and HV components, respectively. On the contrary, at large 
distance from the source, in HH34 as well as in HH1, n e 
increases with velocity. The average values of the electron 
density along HH1 are n e ~ 9.8 x 10 3 cirT 3 and n e ~ 9 x 10 3 
cirT 3 for the blue- and red-shifted components, respectively. 

- The mass flux is mainly carried by the high velocity gas in 
both jets. We derive a lower limit on the mass flux of 3 - 8 x 
10 -8 M yr" 1 along the HH34 and HH1 jets. 

- Comparing the observed [Ti n] /[Fe n] ratio with the theore- 
tical one, we derived a gas-phase Fe/Ti abundance ratio 2-3 
times less than solar. This seems to indicate that the release 
of gas-phase elements from dust grains likely follows a se- 
lective pattern in which Ti -bearing condensates are more eas- 
ily destroyed than the Fe ones. 

- By a comparison of our PV diagrams and electron densi- 
ties with those provided by models for MHD jet launching, 
we suggest that the kinematical features observed close to 
the source in our spectra can be, qualitatively, reproduced by 
both disc-wind and X-wind models, although none of them is 
able to explain the persistency of the LVC at large distances 
(up to 1000 AU) from the launching region. Moreover, none 
of the excitation mechanisms proposed in these models can 
explain the dependence of electron density with velocity that 
we measure in HH34. Alternatively, the LVC of the [Fe n] 
lines that we observe in HH34 could represent dense gas en- 
trained by the high velocity collimated jet. 

- In the H2 line, only the LVC is observed down to the central 
source, while the HVC is detected only up to 2" from the 
star. We suggest that the low velocity molecular gas could 
be excited by oblique shocks occurring along the wall of a 
cavity created by the interaction of a wide angle wind or by 
the external layers of a disc-wind. On the other hand, the 
high velocity component of the H2 emission could be due to 
reformation of H2 molecules along the jet, since the dust in 
the jet is probably not fully destroyed. 
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